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Abstract. We compare performances of ground-based single-mode and multimode (speckle) interferometers in the presence 
of partial Adaptive Optics correction of atmospheric turbulence. It is first shown that for compact sources (i.e. sources smaller 
than the Airy disk of a single telescope) not entirely resolved by the interferometer, the remarkable property of spatial filtering 
of single-mode waveguides coupled with AO correction significantly reduces the speckle noise which arises from residual 
wavefront corrugations. Focusing on those sources, and in the light of the AMBER experiment (the near infrared instrument 
of the VLTI), we show that single-mode interferometry produces a better Signal-to-Noise Ratio on the visibility than speckle 
interferometry. This is true for bright sources (K < 5), and in any case as soon as Strehl ratio of 0.2 is achieved. Finally, 
the fiber estimator is much more robust - by two orders of magnitude - than the speckle estimator with respect to Strehl 
ratio variations during the calibration procedure. The present analysis theoretically explains why interferometry with fibers can 
produce visibility measurements with a very high precision, 1% or less. 
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1. Introduction 

The great interest of using spatial filtering properties 
of optical waveguides in astronomical interferometers has 
been proven in the past years (Coude du Foresto et al. 1997, 
Berger et al. 2001). As a consequence, integrated optics and 
fibers are more and more introduced in the design of 
present and future interferometers to caiTy the signal to 
the detector. Furthermore, practical and theoretical studies 
(Haguenauer 2001, Guy on 2002, Mege et al. 2003) have been 
undertaken to investigate the physical and optical properties in 
waveguided interferometers. 

The present work aims at comparing the sensitivity and ro- 
bustness of single-mode and multimode (speckle) interferom- 
etry. In Section 2, we recall the basic concepts of the fibered 
interferometric equation and the modal visibility. We derive, 
in Section 3, the formal expression of the Signal to Noise 
Ratio (SNR) of the modal visibility which takes into account 
photon, detector and atmospheric noise. In Section 4, we pro- 
pose an analytical approach to estimate the profile of the vis- 
ibility SNR as a function of the magnitude, from partially 
Adaptive Optics (AO) corrected interferograms. We also de- 
rive the performance of single-mode interferometry applied 
to the AMBER experiment (the near infrared instrument of 
the VLTI), in the case of single Gaussian sources. Finally, in 
Section 5, we compute the performances of the multispeckle 



method (Roddier & Lena 1984, Mourard et al. 1994) currently 
used to estimate visibility from non-fibered interferometers and 
we compare the performances and the robustness of single- 
mode and speckle interferometry. 

2. The modal visibility 

Mege (2002) theoretically described how the light is caiTied 
and processed through an interferometer with optical waveg- 
uides. Specifically, he highlighted the coupling phenomenon 
between the incoming wavefront and the fiber, and analyzed the 
characteristics of the interferogram recorded on the detector. 
The principal results concerning the interferometric equation 
of a N telescopes {N,ei) fibered interferometer and the resulting 
modal visibility can be summarized as following: 
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where /(/) is the Fourier transform of the interferogram at the 
spatial frequency /, Vi, is the visibility of the source and Ki 
is the number of photoevents from the telescope that would 
be detected in absence of fibers. //,(./) and - fij) are the 
Fourier transforms of the so-called normalized caiTying waves 
centered at respectively the frequencies / = and / - fjj 
(Mege et al. 2001). Their shape, hence the shape of the inter- 
ferogram, is entirely determined by the geometry of the fibers. 
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p,(y*) and p/yCV*) are respectively the low- (LF) and high- 
frequency (HF) instantaneous coupling coefficients. They give 
the fraction of flux entering in the fibers respectively for the 
photometric and the interferometric channels. Their theoretical 
expressions are: 



Pij(y*) = Po 



V/W)d//r^(/)d/ 



(2) 
(3) 



where V and T'^ are the (partially AO corrected) modal 
transfer functions resulting respectively from the nor- 
malized auto-correlation and cross-correlation of the en- 
trance aberration-corrupted pupil weighted by the fiber 
single mode (Roddier 1988, Mege et al. 2001). T'^ is the 
turbulence-free modal transfer function and po is the op- 
timum coupling efficiency fixed by the fiber core design 
(Shaklan & Roddier 1988). Equations 2 and 3 generalize the 
turbulent coupling efficiency to any kind of sources, i.e. not 
only when it is unresolved by a single telescope (as first no- 
ticed by Dyer & Christensen (1999) in the non-turbulent case). 
Note that, for a point source, the low frequency coupling 
coefficient is proportional to the Strehl ratio S, p - poS 
(Coude du Foresto et al. 2000), and the high-frequency cou- 
pling coefficient follows the simple relationship |py p = pipj. 

From Eq.'s 2 and 3, we can deduce the expression of the 
instantaneous modal visibility Vfj at the spatial frequency fij. 
It is defined as the ratio between the coherent energy (high fre- 
quency) and the incoherent one (low frequency); 
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Perfect equality between the instantaneous modal visibility and 
the object visibility exists only for point sources (V^. = = 
1). In the general case, however, the instantaneous modal visi- 
bility does not match that of the object, especially if the source 
is extended. 

In terms of measurements, the estimator of the modal vis- 
ibility y? can be computed as the ratio between the interfer- 
ogram power spectrum at the frequency fif \I^{fij)\ = (1 - 
T)^\pij(y*)?' KiK j , and the photometric fluxes: ki - Tpi(Vi,)Ki, 
where t is the fraction of light selected for photometry at the 
output of the beam-splitter. Assuming that the latter are esti- 
mated independently through dedicated outputs of the optical 
waveguide (so-called photometric outputs), it holds: 
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In the next Section, we derive a formal expression of the 
relative error - the inverse of the SNR - of the modal visibility. 



also Chelh 1989, Berio et al. 2001). Within this framework, die 
signal is assumed to be corrupted by three difl'erent types of 
noise: (i) the photon noise, with (1 - t)K and tK the total 
number of detected photoevents in the interferometric and the 
photometric channels respectively; (ii) the additive Gaussian 
noise of global variance cr^ which arises from the detector and 
from thermal emission; (iii) the atmospheric noise which re- 
sults from the coupling efficiency variations due to the turbu- 
lence. 

The noise contributions are derived in Appendix A and B. 
The square of the relative error of the modal visibility can be 
described as the sum of three terms: 

= eliK,p) + &\{K, a^,p) + ej (p), (6) 

y2. 

y 

where Sp(K,p), Ba(K, cr^,p) and Ssip), are the relative errors 
due to photon, additive and atmospheric noise, respectively. It 
holds: 
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where X denotes the expected value of the random quantity X, 
and Cov{X, Y] the covariance between X and Y. pif is the aver- 
age over all the telescopes of the LF couphng coefficient, and 
Npix is the number of pixels per interferogram (see Appendix 
B for more details). 

For bright sources, with the exception of point sources for 
which Ipyp = pipj (and hence £s(p) = 0), the dominant 
noise is the atmospheric noise, which results from the classi- 
cal speckle noise filtered by the fiber. From now on, we refer 
to that noise as the modal speckle noise, and the correspond- 
ing SNR, (£s(p)"'), as the modal speckle SNR. This depends 
on the variation of the LF and HF coupling coefficients and 
hence, on the strength of the turbulence. As in the speckle case, 
it does not depend on the source brightness and corresponds to 
the maximum achievable SNR per interferogram. Nevertheless, 
as we will show, it depends on the source size. 



3. Relative error of the modal visibility 

The noise calculations are based on the spatially continuous 
model of photodetection introduced by Goodman 1985, (see 



4. Performances of single-mode interferometry 

In this Section, we develop a simple model to estimate the 
modal speckle noise from partially AO corrected interfero- 
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Fig. 1. Modal SNR per interferogram as a function of the magnitude of the source in the K band. From left to right, top to bottom, the object 
size increases: point source, 1/50, 1/20, 1/10 of Airy disk. The considered baseline is Bx = By = 100m which mimics an average baseline of 
the VLTI array. This corresponds to visibility of respectively 1, 0.882, 0.457, 0.04. Results are given for different Strehl ratios: S = 0.011 (solid 
line, no correction), 5 = 0.1 (dashed line), S = 0.5 (dash-dotted Une) and S = 0.9 (dotted line). 



grams, and we derive the performances of single-mode inter- 
ferometry applied to AMBER, the near-infrared instrument of 
the Very Large Telescope Interferometer (VLTI). 



4.1. The modal speckle noise 

The derivation of the modal speckle noise is fully detailed in 
Appendix C: here we recall the outlines of our approach. We 
assume that: (i) the distance between the telescopes (baseline) 
is larger than the outer scale of the turbulence, i.e. that the com- 
plex ampUtudes over pupils i and j are uncorrelated; (ii) the 
atmospheric phase has Gaussian statistics; (iii) the associated 
structure function is (spatially) stationary. 

At this stage, the modal speckle SNR is described by in- 
tegrals of dimensions up to sixteen (see Table C.2, from Eq 
C.ll to Eq C.15), which depend on the pupil, the object vis- 
ibiUty and the structure functions, respectively. Here we esti- 
mate these integrals with a simplified analytical approach deal- 
ing with Gaussian functions only. We assume that the object 
brightness distribution 0,(a) and the pupil function have re- 
spectively the form: 



0,ia) = exp — T- 



(10) 



and, 
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where D is the diameter of the telescope and o-p is defined in 
such a way that the integral of the Gaussian pupil, i.e. its sur- 
face, is equal to the integral of a circular pupil of diameter D. 
Furthermore, in order to take into account the partial correction 
by adaptive optics, we approximate the structure function 2)(/) 
by: 

1"-' h + il-h)!B<!,{f) (12) 



exp 



where h € [0, 1] defines the level of correction and So(/) is 
the transfer function of the turbulent atmosphere (Conan 1994) 
that we will assume Gaussian too: 



Sfl) = exp 



(13) 



ro being the Fried parameter. Note that, from the previous equa- 
tions, the Strehl ratio S is given by: 



S = h + {l-h)- 



(14) 



At this point, the integrals C. 13, C. 14 and C.15 caimot be com- 
puted formally yet. To do that, we perform a limited expansion 
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Fig. 2. Instantaneous limiting magnitude of single-mode interferome- 
try technique, defined such that the modal SNR per interferogram is 
equal to unity. From top to bottom we plot limiting magnitudes for ob- 
ject sizes of respectively: point source, 1/50, 1/20, and 1/10 of Airy 
disk. Note that Eq. 15 reproduces fairly well the curves for object sizes 
equal or smaller than 1/20 of Airy disk. 



of the expressions under the integrals with respect to the vari- 
able e = (1 - K)lh. We end up with a series development in 
e" in which the coefficients are integrals of Gaussian functions 
products only. These coefficients are then computed with the 
MAPLE software. However, the series development converges 
only for e e [0, 1[ or h e]0.5, 1], i.e. for good to perfect AO 
corrections (see Eq. 14). In order to estimate the modal speckle 
SNR for average to low AO correction, we compute it in the 
pure turbulent case (h - 0) and we extrapolate from /z = to 
h = 0.5. The case /i = is computed separately from stan- 
dard hypothesis assuming that the complex amplitude of the 
pure turbulent wavefront follows circular Gaussian statistics 
(Roddier 1988). We finally obtain an expression of the modal 
speckle SNR as a function of the major parameters of the ob- 
servation: the source size, the baseline, the turbulence strength 
D/ro, and the level of AO correction, i.e. the Strehl ratio S. 
Strictly speaking, the structure function is not stationary as the 
error of the AO corrected wavefront increases from the center 
to the edge of the telescope pupil. Also, its shape is not ex- 
actly described by the simplified Eq. 12. Nevertheless, we ex- 
pect the modal speckle SNR resulting from our model, deaUng 
with Gaussian functions only, to have the right order of magni- 
tude and the correct functional dependencies. 

4.2. Performances 

In this Section we compute the SNR profiles of the modal vis- 
ibiUty (per interferogram) for faint to bright compact sources 
(i.e. smaller than the Airy disk of a single telescope). We con- 
sider observations under average seeing conditions (ro - 1 .6m) 
with 2 Unit Telescopes (D = 8m) at the baseUnes Bx - By = 
100m and the AMBER recombiner in the K band {2.2nm). We 
assume that the interferogram is dispersed along the columns 
of a bi-dimensional detector and that each spectral channel is 
sampled with 6 pixels to ensure low and high frequency peaks 



separation. For quantitative calculations, we choose the spe- 
cific instrumental parameters of AMBER (Malbet et al. 2000, 
Petrov et al. 2000) together with a spectral resolution of 35, an 
integration time of 30ms per interferogram, a detector read- 
out noise of 15e~/pixel and a transmission coefficient r = 0.5. 
Note that, in those conditions, thermal noise is negligible. We 
also assumed an optimized instrumental coupling efficiency 
po = 0.8 (Shaklan & Roddier 1988). Fig. 1 shows for different 
Strehl ratios, the modal visibility SNR as a function of the mag- 
nitude for 4 object sizes, point source, 1 /50, 1 /20 and 1/ 10 of 
Airy disk, with respective visibiUties 1, 0.882, 0.457 and 0.04 
at the baselines Bx - By - 100m. 

We can clearly see that the saturation regime, where the 
modal speckle noise dominates, is significant only in the ab- 
sence of AO correction. As soon as the image is partially AO 
corrected, even at small Strehl ratios {S > 0.1), the satura- 
tion regime is rejected towards negative magnitudes. It is re- 
placed by an extended 'photon noise' regime, which depends 
on the total flux weighted by the statistics of the coupling co- 
efficients, that can be interpreted as transmission coefficients. 
These transmission coefficients decrease with the Strehl ratio 
and also when the size of the object increases, therefore lower- 
ing the modal visibility SNR. At low fluxes, the 'detector noise' 
regime, marked by the break in the SNR slope, takes over. 

We define the instantaneous limiting magnitude as the mag- 
nitude for which the SNR per interferogram is equal to unity. It 
occurs at very low fluxes, where the additive noise is dominant. 
From Eq. 8, and after some simplifications which are valid for 
partially resolved objects (typically ag > 1/20 of Airy disk 
at the considered baseUne), the corresponding limiting flux is 
given by: 



Klim — 



xCrNtel 



il-T)VijPoS 



(15) 



Figure 2 shows the limiting magnitude as a function of the 
Strehl ratio for four object sizes: point source, 1 /50, 1/20 and 
1/10. Without AO correction, the limiting magnitude is small, 
between 5 and 6. As soon as the image is AO corrected and the 
object partially resolved, the limiting magnitude significantly 
increases, reaching about 10 for a Strehl ratio of 0.5. However, 
for largely resolved objects (^ 1/10 of Airy disk), the cou- 
pling efficiency becomes so low that performances of fibered 
interferometers in terms of SNR and limiting magnitude are 
severely degraded. 

Note, however, that calculations (their detailed description 
is beyond the scope of this paper) show that the saturation 
regime can span higher ranges of magnitude in cases where 
the compact Gaussian source is surrounded by an extended dif- 
fuse matter such as a disk or a dust shell. Indeed, it can be 
shown that, depending on the fraction of the flux in the ex- 
tended structure, the modal speckle SNR can decrease by an 
order of magnitude, or more. This is due to the strong fluctua- 
tions of the LF coupling coefficient of the extended component. 
Such behavior is specific to single-mode interferometry, and it 
demonstrates that fibers are more efficient when they deal with 
compact sources. On the contrary, we emphasize that in classi- 
cal (multimode) interferometry, a diffuse extended component 
has no effect on the SNR. 
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Fig. 3. Multispeckle SNR per interferogram as a function of the magnitude of the source in the K band. As Fig. 1, which those curves must 
be compared to, the object size is varying from left to right, top to bottom: point source, 1/50, 1/20, 1/10 of Airy disk, for a baseline 
Bx = By = 100m. Results are given for different Strehl ratios: S = 0.011 (solid Une, no correction), S = 0.1 (dashed line), S = 0.5 
(dash-dotted line) and 5 = 0.9 (dotted Une). 



5. Single-Mode versus multimode interferometry 

In this Section, we compute the performances of speckle inter- 
ferometry and we compare their performances and their robust- 
ness to those with single-mode interferometry. 

5.1. Performances of speckle interferometry 

In absence of waveguides, visibility estimators can be de- 
fined by using speckle techniques, following Labeyrie's 
method (Labeyrie 1970, Sibille et al. 1979) for single dish 
observations (see Appendix D). The classical estimator 
(Roddier & Lena 1984) consists in taking the ratio of the in- 
tegral of the high frequency spectral density by the integral 
of the low frequency one. Taking the integral of the high fre- 
quency peak is essential to perform a consistent comparison 
with the fiber case, since fibered interferometry induces an av- 
erage (more precisely a convolution) of the visibility over the 
high frequency (see Eq.3). Moreover, to insure a thoroughly 
consistent comparison, we assume that the photometric fluxes 
are measured simultaneously with the interferograms (instead 
of taking the integral of the low frequency spectral density) and 
we define the estimator of visibiUty as the ratio between the in- 
tegrated high frequency peak of the spectral density and the 



photometric fluxes. It may be written as: 

_ J |/(/)l'd/ 

This estimator needs to be calibrated by a point source. It 
has basically the same performance than the classical speckle 
estimator, but not the same robustness, as we wiU see later. 
In practice, we replace the integral by a discrete sum with a 
regular spacing A/ - ^. Since we consider partial correc- 
tion by Adaptive Optics, the points involved in the discrete 
sum are not statistically independent. Their correlations are 
taken into account in the signal to noise ratio calculations, de- 
tailed in Appendix D. These calculations require the knowl- 
edge of the first and second order statistics of the speckle trans- 
fer function. Formal expressions of those moments have been 
derived following a procedure similar to the one described in 
Section 4.1. For numerical appUcations, we adopted the same 
parameters than for single-mode interferometry, with the ex- 
ception of the number of pixels required to correctly sample 
the interferograms, which is = 6(D/ro) for 2 telescopes 
(ChelU & Mariotti 1986). Fig. 3 shows, for different Strehl ra- 
tios, the speckle visibility SNR as a function of the magnitude, 
for the four object sizes previously considered. The striking dif- 
ference with single-mode interferometry is the almost absence 
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Fig. 4. Ratio between fibered and multispeckle visibility SNR as a function of the magnitude for different AO correction levels, from left to 
right, top to bottom: .S = 0.011 (no correction), S = 0.05, S = OA,S = 0.2, S = 0.5, S = 0.7, S = 0.9 and S = 1. Ratios are displayed for 
different source sizes in fraction of Airy disk: point source (solid line), 1/50 (dashed line), 1/20 (dash-dotted line), 1/10 (dotted line). 



of the 'photon noise' regime, even with AO correction. Instead, 
the saturation regime is reached in the entire range of magni- 
tudes until the 'detector noise' regime takes over. Note, further- 
more, that the speckle noise is barely dependent on the source 
size and does not cancel out for a point source, contrarily to the 
modal speckle noise. 



5.2. Comparison of performances and robustness 

Fig. 4 shows the ratio between single-mode and speckle SNR as 
a function of the magnitude, for different compact source sizes 
and Strehl ratios. Clearly, with the exception of bright sources 
(K < 5), speckle interferometry without AO correction is su- 
perior to single-mode interferometry. However, as the degree 
of AO correction increases, the situation very rapidly evolves, 
and, starting from a Strehl ratio of 0.2, single-mode interferom- 
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Fig. 5. Comparison of flbered and speckle interferometry limiting 
magnitudes. Limiting magnitude is still defined such that SNR - 1. 
Solid lines and dashed lines show fiber and speckle cases respectively. 
From top to bottom, the object size is: point source, 1/50, 1/20, and 
1/10 of Airy disk. 

etry always reaches higher SNR than speckle interferometry. 
Such efficiency of fibered interferometry is due to two major 
aspects: the remarkable property of spatial filtering of the fiber 
as soon as enough coherent energy (typically 20%) lies in the 
Airy disk, together with the possibility of sampling the signal 
on few pixels. On the contrary, multimode interferometry re- 
quires at least good AO correction (S > 0.5) to significantly 
reduce the speckle noise. Moreover, fringe sampling is seeing 
dependent and requires a much larger number of pixels, spe- 
cially for telescopes with large apertures such as the VLTI. Fig. 
5 compares the limiting magnitudes of both methods. Without 
AO correction, the speckle limiting magnitude is between 8 and 
9, well above the single-mode one. With AO correction and 
Strehl ratios larger than 0.2, the limiting magnitudes are basi- 
cally the same. 

To evaluate the robustness of each method, we investigate 
the stability of the measured visibility versus the Strehl ratio. 
Fig. 6 shows the visibility variation as a function of the Strehl 
ratio, normalized to the visibility at .S = 0.7. The two upper 
curves correspond to the speckle estimator studied above and 
to the classical speckle estimator (ratio of the integrated high 
and low frequency peaks of the spectral density), respectively. 
The classical speckle estimator is much more robust, but, in 
both cases, even a small Strehl ratio variation (0.2) can pro- 
duce visibility variations up to 10%. The 3 lower curves corre- 
spond to the modal visibility estimator for source sizes 1/50, 
1 /20 and 1 / 10 of Airy disk (point source is irrelevant since it is 
theoretically independent of the turbulence). The robustness of 
the modal estimator depends on the source size, while the mul- 
tispeckle estimator does not. In any case, however, the modal 
estimator is clearly more robust than the speckle one, by more 
than 2 orders of magnitude. The modal visibility is stable at a 
level less than 1% over all the range of possible Strehl ratios, 
from to 1. This last property is interesting, not only to per- 
form high precision measurements, but also for the selection 
of reference sources. Indeed, it would suggest that the use of a 
reference source having a large magnitude difference with the 
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Fig. 6. Robustness of single-mode and multimode estimators as a 
function of the Strehl ratio. Comparison is given in terms of relative 
error on the visibility AV/V. Solid lines deal with fiber estimator for 
respectively 1/10, 1 /20 and 1 / 100 of Airy disk, top to bottom. Dashed 
lines gives the robustness of two speckle estimators. Top: so-called es- 
timator I which takes into account the estimator defined in this paper 
(Eq. 16). Bottom: estimator U describes the "classical" speckle es- 
timator (see text). Note that both speckle estimators are object size 
independent. 

studied object, and hence, where the level of AO correction is 
different, does not affect the precision of the measurement. 

Note, however, that this comparison focuses on compact 
sources. For a central source surrounded by a larger diffuse 
component, the fluctuations of the LF coupling coefficient due 
to the extended structure severely reduces the filtering prop- 
erties of the single-mode fibers. As a consequence, a situation 
similar to the multispeckle case occurs with a saturation regime 
spanning a large range of magnitudes, and where the larger the 
fraction of flux in the extended component, the larger the range. 
The presence of an extended structure also causes a decrease of 
the robustness, which, however, remains better than in the case 
of speckle estimators, by one order of magnitude. 

6. Summary 

In this paper, we have developed a formalism that can predict 
theoretical SNR on visibilities, when measurements are par- 
tially AO corrected and are corrupted by photon, additive (de- 
tector, thermal) and residual atmospheric noise. This formalism 
has been applied to single-mode and multimode (speckle) in- 
terferometers. We have assumed that: (i) the wavefronts over 
two distinct telescopes are uncorrelated; (ii) the atmospheric 
phase has Gaussian statistics, and the associated structure func- 
tion is (spatially) stationary; (iii) the pupil function, the object 
brightness distribution, and the transfer function of the turbu- 
lent atmosphere are Gaussian. In the case of single-mode in- 
terferometry and for compact sources (i.e. sources smaller than 
the Airy disk of a single telescope) not entirely resolved by 
the interferometer, we show that, in the presence of AO cor- 
rection, the remarkable filtering properties of fibers reject the 
SNR saturation regime due to speckle noise towards negative 
magnitudes. Instead, the modal visibility SNR is dominated by 
the 'photon noise' regime followed by a break close to the lim- 
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iting magnitude due to 'detector noise' . For AMBER, the near 
infrared experiment of the VLTI, we found that the limiting 
magnitude is about K = 10 for a Strehl ratio of 0.5, with two 
8m Unit Telescopes (typical atmospheric and instrumental pa- 
rameters of the VLH/AMBER instrument were discussed in 
Section 4.2). In the case of speckle interferometry, the visibil- 
ity SNR is characterized by the almost absence of the 'pho- 
ton noise' regime, replaced by the saturation regime due to 
speckle noise. With the exception of bright sources (K < 5), 
without AO corrections speckle interferometry reaches higher 
SNR than single-mode interferometry, leading to higher limit- 
ing magnitudes. However, the situation is different in presence 
of AO correction, even small (S - 0.2): the hmiting magni- 
tude of both methods becomes the same, but the single-mode 
reaches much higher SNR than speckle interferometry. This 
last property, together with the insensitivity of the modal visi- 
bility to the Strehl ratio variations, explains why interferometry 
with fibers can produce visibility measurement with a very high 
precision of 1% or less, on compact sources. 

Appendix A: Moments of the spectral density of 
an interferogram 

A.1. Formalism 

In order to compute the moments of the spectral density, we 
use the spatially continuous model of photodetection process 
of Goodman (1985) where the detected signal takes the form: 

K 

dix,y) = ^6(x- x„,y-y„) (A.l) 

n=l 

and its Fourier transform: 

K 

D(f,, fy) = Yj c-^''^""^^''"^ (A.2) 

The position (x„, y„) as well as the number of photoevents K are 
considered as independent random processes with probability 
density functions proportional to the intensity Iix,y). This for- 
malism has been deeply studied (Goodman 1985, Chelli 1989) 
and the results already proven will not be re-demonstrated 
here. This appendix especially focuses on the original approach 
where the signal is corrupted by an additive noise e, following 
a centered Gaussian random process, and of Fourier transform 
T. Using a mono-dimensional writing, but without loss of gen- 
erality, the corrupted signal S (/) has now the expression: 

K 

S(f) = D(f) + 7(f) = Yj e"'"^'" + "^Z) (A.3) 

n=l 

We compute the spectral density statistics with respect to the 
additive Gaussian noise, the photon noise and the turbulent at- 
mosphere, respectively. The estimator of the spectral density 
is: 

Qif) = \S{f)f (A.4) 



A.2. Expected value 

We first take the expected value with respect to the additive 
noise from Eq. A.3. It becomes: 

E.{e(/)} = |5(/)|2 + EJ[^2j ^ |5(^)|2 + (^ 5) 

Then computing |Z)|^ from Eq. A.2: 

K K 
n=l 1=1 

Goodman formalism finally leads to a total expected value of: 

E{e(/)} =< >^ + <K>^ +N(r' (A.7) 

where denotes the expected value relatively to the turbu- 
lent atmosphere random process. 

A.3. Variance 

The variance is defined by: 

(rHQ}^E{Q^}-E{Qf (A.8) 

Computing E{2^) with respect to the additive noise, we have: 

E{e^) = E{\D\'^} + 4NcrM\Df] + 2N^cr'^ (A.9) 

from which Goodman formalism in the pure photon noise case 
can be appUed. Finally the variance of the spectral density may 

be written as: 

cr^iQif)} =< ^R/)i' - < K^mf >l + 

4 < K^Viff >^-2<K >^< K^lKff >^ + 

4 < K^\i(.f)f >^+2<K^>^-<K >l +3N(T^ + 

N^cr'^ + 2Na^ <K>^ +2Ncr^ < K^ViDf ><l> (A.IO) 

A.4. Covariance 

We can also compute the covariance of the spectral density 
which is defined as: 

Coviei.gi} = E{QiQ2]-E{Qi}E{Q2] (A.ll) 

It becomes: 

cov{e(/i), e(/2)} =< 'tiiifitmt >^ - 

< K^mt >^< '^mt >0 + 

< tiifS'ifiWifi - fi) >0 + 

<'tf{f,)i{f2)i{fl-fl) >0 + 
<Z^(/l)T(/2)?(/l+/2) >0 + 

< :^^*(/i)?(/2)T(/i + /2) >^ + 

< K'mt >^-<K >^< 'tiiifit >^ + 

< K'vifit >^-<K >^< 'tmt >^ + 
2 < 'tmf >0 +2 < r'R/2)i' >0 + 

< K'm - fit >0 + < K'm + fit >0 + 

<K^>^-<K>l + <K>^ +3Ncr'^ (A.12) 
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Appendix B: Estimator of the visibility for fibered 
interferometers 

B.1. The interferometric equation 

Mege (2002) has shown that the fibered interferometric equa- 
tion could be written as follows: 

/(/) = J]KiPiHtif) + 

i 

where p, and pij are respectively the low and high frequency 
coupling efficiencies and Ki is the number of detected photo- 
events in the absence of a fiber on the i''' telescope. Moreover 
we define ^', as the number of photoevents at the output of the 
fiber. We have the relationship ki = piKi. For simplicity pur- 
poses, we also introduce a "global" low-frequency coupling co- 
eflicient p;/ such as: pi/K = Yjj PjKj where K is the total num- 
ber of photoevents, i.e. K = XjKj. Supposing that V„Pi = p, 
wehavep;y =p,- =p. 

B.2. Estimator of tiie visibility/Error on ttie modal 
visibility 

A classical estimator of the modal visibility consists in divid- 
ing the HF spectral density of the interferogram by the photo- 
metric fluxes. We assume in the following that a fraction t of 
the light has been selected for photometry analysis and that the 
rest of the light (1 -t) belongs to the interferogram. For sake of 
simplicity, we define and as the total number of photo- 
events concerning respectively photometric and interferometric 
channels, i.e K"^' = tK and = (1 - t)K. It comes: 



(B.2) 



Using Papoulis (1984) second order approximation, we can de- 
rive the square variance relative error of the modal visibiUty: 



—2 2 



coY{mj)\\kfk^} 



-2 



E{mj)mfk^ 



(B.3) 



To derive those moments, we use the Goodman formalism de- 
scribed in Appendix A. In the case of fibered interferome- 
ters, the LF and HF of the Fourier-transformed interferogram 
(Hiif) and Hijif)) are fixed by the geometry of the fibers, and 
statistics with respect to the turbulent atmosphere only appear 
in the coupling coefficients. For sake of simplicity we assume 
those coupling coefficients to be uncorrected between two dif- 
ferent baselines, although it does not change the eventual con- 
clusions. Assuming also that the telescope transmissions are all 
equal, i.e. Ki = K/Ntd, it leads to the following expressions: 



\pij\2—+pKi + NcT' 



(B.4) 



Tel 



[4pif\pij\2-2m^ — 



4boP— +2Arcr2boPT^ + 



<r\k^k^,) 



tel 



tel 



tel 



li^P^ r 

Cov{|/(i;-^)p,^f^f} = -——\\p,^\^p,p^ - \p]0' 



(B.5) 
(B.6) 

(B.7) 



And the square relative error on the modal visibihty can be 
expressed as the sum of three contributions: 



o-2{y2.} 

= &l{K,p) + &liK, ali,p) + &l(p) 



(B.8) 



y2 
ij 



where: 



N,A'^p,f\Pif-2p\pif) 



(1-T)h/ 



Ay 
K 

I 



P 



(B.9) 



(l-ry\Pif K 
is the photon noise square relative error: 



(I-T)V/ ^ 0--TY\pij\'K 

INo-^p 



^2 —3 



(B.IO) 



takes into account the additive noise cr^ and 

p2_ ^'^liPj C0V{p2.,p;p,.} 
t>c — r + ;; Z- 



\Pij\- 



(B.I I) 



arises from the couphng fluctuations inducing a so called 
"modal-speckle" square relative error Concerning the latter 
contribution, it is interesting to notice that in the case of a point 
source, we have \pij\^ - pipj. The outcome is that the "modal- 
speckle" noise contribution for a point source is zero. 



Appendix C: Coupling coefficients statistics 

In order to calculate the square relative error of the modal vis- 
ibility given in Appendix B, we first compute first and second 
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Table C.l. LF and HF coupling coefficient first and second order statistics. Rigorous expressions. 



pi = ^ J Vl(u)Pia)Pia + u)<'i'iia)%ia + u)>dadu (C.l) 

Ia/ = ^ J K(« - fij)V^(v - fij)P(a)P(a + u)P(fi)P(fi + v) < ^,(a)'i'l(J3)^}(a + u)'i'j(J3 + v) > dadM^Sdv (C.2) 

E{|pi/} = 

J V*(fij - aiWlifij - a2)VMi - a^Wlifij - a,)P{a)P{a + ai)P(fi)P(fi + a2}P(y}P(y + a3)P{d}P(S + a,) 

< %{a)'¥;(fi)'¥i{y)'¥'{d)'i"j{a + ai)'Pj(j8 + aiWjiy + a^Wjid + 04) > daida2da3da4dadfidydS (C.3) 

J K(ai}V:{a2)V:{a3)V:{a4)P(a)P(a + ai)P(J3)P(J3 + a2)P{y)P{y + ai)P{6)P{8 + a^) 

< 1',(ar)*P*(a + ai)1'i(j8)*'Pi(j8 + 02)'^' j{yyv*{y + a^fV j{6yv){6 + 04) > daida2da3da4dadfidyd6 (C.4) 

ji j V*(fij - aiWtifij - a2)V:{a3)V:{a4}P(a)P(a + ai)P(fi)P(fi + a2)P{y)P{y + a3)P{d)P(S + 04) 

< %(a)'¥;(fi)'¥i{y)'¥i{y + a3)1'*(ar + aO'i'jifi + a2Wj(S)'i"j{S + 04) > daida2da3da4dadfidydS (C.5) 



order moments of the LF and HF coupling coefficients which 
are defined as following: 



Pi(V*) = po 



P,;(^*) = Po 



iV**T')f=o 
^jT',(f)dfjT^(f)df 



that (Conan 1994): 

D(u, v) =< ij/i(u)ij/*(v) >^=< j{u)\l/){y) >^ 



(C.IO) 



(C.6) 
(C.7) 



Finally we suppose that these structure functions are stationary, 
i.e. 2)(m, v) = 2)(m - v). After such formal derivations we obtain 
expressions of the coupling coefiicients statistics as summa- 
rized in Tables C.2. Readers may note that in the pure turbulent 
case, we assume in addition that the complex amplitude of the 
wavefront follows circular Gaussian statistics (Roddier 1988), 
hence slightly changing the expressions given in table C.2 

Appendix D: Estimator of the visibility for 
multispeclde interferometry 

D. 1. The inteiierometric equation 

For non fibered interferometers, the convolution between the 

object and the interferometer transfer function stands for any 
spatial frequency (Tallon & Tallon-Bosc 1992). We may ex- 
press: 



where T' and T'^ are the (partially AO corrected) modal trans- 
fer function resulting respectively from the normalized auto- 
correlation and cross-correlation of the entrance aberration- 
corrupted pupil weighted by the fiber single mode, i.e: 

T'iu) = J j Piia)Piia + u)'i'iia)Xia + u)da (C.8) 

Introducing Eq.'s C.8, C.9 in Eq.'s C.6 and C.7 and developing 
the expressions of first order (pi, |p,jp), and second order (f^.p. 

cj^.p., Cov{|p,/,p/Pj)) coupling coeflicient statistics, leads re- j^^jr^ ^ y ^.-^y-, + V V JWiKjiif - fij) (D.l) 

spectively to second, fourth and eighth order moments of the ^ it 

complex amplitude of the wavefronts ^iiu) and ^P/m). Such 
rigorous expressions are written in Table C.l. 

At this point, those expressions are not yet formally com- 
putable. Hence we perform simplifications of the equations. 
We first assume that the wavefronts are uncorrelated, i.e. < 
'¥i(u)'¥jiu) >^=< »Pi(M) >^< *Py(M) >^. Then we use Korff's 
(1973) derivation of the moments of the complex amplitude of 
the wavefronts to introduce in the equations, linear combina- 
tions of the structure function (Diu, v)) at different spatial fre- 
quencies. We recall that the structure function is defined such 



with 



;■(/) = y(/).5(/) 



(D.2) 



where V(f) is the object visibility and S (/) is the normalized 
transfer function of the interferometer which consists in the au- 
tocorrelation of the pupil function weighted by the remaining 
atmospheric phase. We have the following relationship 

s{f) = J]s '(/) + ^ 2 Z ^ '''^^ 
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Table C.2. LF and HF coupling coefficient first and second order statistics, simplified expressions after assuming: (i) no correlation between 
wavefronts; (ii) KorfF's derivation of the moments of the complex amplitude of the wavefronts; (iii) stationarity of the structure function 



pl = ^ J v:(M)P(a)P(ar + M)e-J^"'dadM (C.ll) 
\pi^ = ^ Jviiu-fijW^iv- fij)P{a)P{a + u)P{p)P{p + v)e- J ^^-'>-^-*^-''''^AaAuApdv (C. 12) 



E{|p,/1 



/ + «i ~ a2)Vl(fij + 03 - a,)V,{f,j + a, - a^WKf.j + «7 - ai)P(ai)P{a2)P{a3)P{ai)P{ai)P{a(,)P{ai)P{ai) 

" p-^lB(fl3-<l5)+0(<I4~a6)+£)(<li-fl7)+£)(a2-<I8)] " 



g-i[B(ai-a3)+B(a5-07)+O(a2-a4)+B(06-<ig)J 



g-i[£)(ai-as)+£)(a2-a6)+B(a3-07)+B(<i4-<is)] 



dfli . . . dfls 



Pi 



J V:(a2 - ai)V:(a4 - UiWlioi - asWlia^ - ai)P{ai)Pia2)Pia3)Pia^)P{a5)P(a(,)P{ai)P{ai) 

g-i[B(<i2-a3)+B(ai-a4)+B(<i6-<i7)+B(a5-a8)]" 



g-^[B(a2-ai)+B(a4-a3)+0(a6-a5)+0("8-"7)l 



E{|p,./pip^} : 



Pi 



g-|[B(a3-ai)H-B(a4-a2)+B(a5-a7)+B(ag-a6)] 



dfli . . . dag 



5-* 



J + ai - a2)V:(/7 + 03 - a^)Vl{at - ai)Vl{ai - an)P{ai)P{a2)P{.ai)P{a^)P{ai)P{at)P{ai)P{a{) 

g-^[B(a3-a5)+B(a6-ai)+B(a7-a2)+B(as-a4)] 



-^[fl(a3-<ii)+fl(a2-a4)+B(<i6-<i5)+fl(a8-a7)] 



g-i[B(a5-ai)H-B(a6-a3)+B(a7-a4)+B(ag-a2)] 



dfli . . . dag 



(C.13) 



(C.14) 



(C.15) 



where and are respectively the auto and cross- 

correlation of single phase corrugated pupil (Roddier 1988), 
i.e.: 

S"{f) = ^ J P{a)P{a + fyi>i{ay¥*{a + f)da (D.4) 
S''{f) = ^ J P{a)P{a + fr¥i{a)'V){a + f)Aa (D.5) 

D.2. Estimator of the visibility/Error on the visibility 

To perform a consistent comparison within the fiber case, we 
define the estimator of the multispeckle visibility as the ratio 
between the integrated high frequency peak and the photomet- 
ric fluxes. We still assume that a fraction t of the light is in- 
jected in the photometric channels. It writes: 



vHfij) ' 



i/(/)rd/ 



' J 



(D.6) 



From a digital point of view, as the size of the instantaneous 
speckle image is ^ wide (perfect correction apart), the integral 
of the previous equation can be replaced by a discrete sum over 
points of the peak that are distributed every ^. We then have 



an equivalent estimator: 



-,A/: 



ro 
A 



We derive the square relative error from Papoulis (1984): 

Z Cov{|/(m)p|/(n)|2} 

^ y'iil I m.n*m 



(D.7) 



n2 



1 1 



(D.8) 



Note that correlations between points have to be taken into ac- 
count since we consider partial correction by Adaptive Optics. 
We then use the formaUsm developed in Appendix A to derive 
D.8. Note that in the multispeckle case, at the contrary of the 
fibered case, the LF and HF peaks of the Fourier-transformed 
interferogram (i.e. 5'(/) and >S'^(/)) depend on the turbulence. 
Moments of |/(/)P are: 



E{|/(/)l'} = < >Vf+Kr + Ncr' 



(D.9) 
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<\s'W ><!,-< \s'W>l vf + 



3Na' + N^cr' 



Cov{|/(m)p|7(n)n = 

Nf ^ 



(D.IO) 



tel 



2 < \S'Kn)\^ >^ Vl+ < \S'{m - n)\^ >^ Vl_„] + 



2— < S'\m)S'\n)S' {m-n)> V„V„V:,_, + 



N: 



tel 



<\S'\m)\'\S'\n)\' 



(D.ll) 



Another classical estimator of the visibility 
(Roddier & Lena 1984, Mourard et al. 1994), known to 
be more robust to the turbulence, consists in taking the ratio 
of the high-frequency peak energy by the low-frequency one 
(coherent energy versus incoherent one): 



mrdf 



(D.12) 
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Assuming that photometric chaimels contain at most 50% of 
the information, we expect the relative error of the visibiUty 
from this estimator to be equal to the previous one within a 
factor of ^/2. 



Appendix E: Speckle transfer function statistics 

In order to calculate the square relative error of the visibility, 
we first need to compute expected values of \S'(f)\^, \S'-'(f)\^ 
and IS'-'if^. As discussed in Appendix D the functions S'if) 
and S'if) are respectively the auto and the cross-correlation 
of the pupil functions weighted by the remaining atmospheric 
wavefront. If we neglect the weighting of the pupU by the sin- 
gle mode of the fiber, we find the very same expression as in 
Eq's C.8, C.9, i.e. S'(u) = r(u) and 5'-'(m) = T'Hu). Deriving 

< \SW >, < \SW >, < \SW > , < \S'W > in- 
troduces once again second, fourth and eighth order statistics 
of the complex amplitude of the wavefronts, that we develop 
following the same procedure as described in Appendix C. We 
finally obtained expressions given in Table E. 1 . 
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< 15 ■(/)!' > = 5-M P(ai)P{ai + a2)P{ai + f)P{a, +a2+f) 



p-[0(fl2)+0(/)] 



-ilffl(fl2+/)+B(a2-/)] 



dflida2 



<\S"{ff> = S-^ J P(ai)P(a2)P(ai+/)/'(a2+/)e-®<'""'"^*daida2 



(E.1) 
(E.2) 



< is'm' >= 

S-^ J PiaOPiai + f)P{a2)P{a2 + f)P(ai + «3)^(ai +03 + f)P{a2 + a^)P{a2 +0^ + /) 

r Q-[Diai-a2)+!Dia3+ai-a2)+2)(n2+a4-ai)+Dia4+a2-ai-a^)] " 



